Universality and universal finite-size scaling functions 
in four-dimensional Ising spin glasses 
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I. INTRODUCTION 

The concept of universality, according to which the 
values of different quantities, such as for example criti- 
cal exponents, do not depend on the microscopic details 
of a model, is well established in the theory of critical 
phenomena of systems without frustration. However, for 
spin glasses,— which have both frustration and disorder, 
the situation has been less clear until recently. Large- 
scale simulations in three space dimensions 2 *^ for dif- 
ferent disorder distributions of the random interactions 
between the spins have shown that the critical exponents, 
the values of different observables at criticality, as well 
as the finite-size scaling functions — all which are nec- 
essary ingredients to show that different models are in 
the same universality class — agree well within error bars. 
The conclusion that universality holds has also been ob- 
tained via other methods such as high-temperature se- 
ries expansions^ where the critical exponent 7 has been 
studied. Studies of dynamical quantities however have 
yielded different critical exponents for different disorder 
distributions^^ although it is unclear up to what level 
it can be expected that "dynamical universality" can be 
compared to universality in thermal equilibrium. 

In two space dimensions things are less clear: Because 
the transition to a spin-glass phase only happens at zero 
temperature, it is believed that systems with discrete and 
continuous coupling distributions behave in a different 
wayJ^ In particular, the ground-state entropy is nonzero 
in the former while it is zero in the latter. Thus one 
might expect that at the zero-temperature critical point 
the critical exponent r\ is different. Recently, however, 
an alternate scenario for universality in two space di- 
mensions has been proposed^ where it is expected that 
as long as the temperature is nonzero two-dimensional 
spin glasses with different disorder distributions belong 
to the same universality class and where the different be- 
havior seen at zero temperature is explained by the pres- 
ence of an additional spurious fixed-point Limitations 
in the simulation techniques and analysis methods have 
so far yielded no conclusive results making this proposal 
controversial J&ll 



In spin glasses it is extremely difficult and numerically 
very costly to determine critical exponents with high pre- 
cision. This is mainly due to the following reason: It is 
difficult to sample the disorder average with good enough 
statistics, especially for large system sizes, because spin 
glasses have very long equilibration times in Monte Carlo 
simulations and as a consequence one has to deal with 
corrections to scaling due to a very limited range of sys- 
tem sizes at hand. In previous studies^ only statistical 
error bars had been considered. Because of limited sys- 
tem sizes and corrections to scaling in three space dimen- 
sions, deviations between the critical parameters beyond 
statistical error bars can be expected and indeed this 
expectation has very recently been confirmed in Ref. |j| 
in a very thorough study where for the first time cor- 
rections to scaling have been studied with good accu- 
racy. While studying higher-dimensional systems than 
three space dimensions might seem paradoxical at first 
because of the aforementioned problems, the proximity 
to the upper critical dimension d UC( j = 6 is advanta- 
geous. High-temperature series expansion studies^ sug- 
gest that corrections to scaling should be falling off fast in 
four-dimensional Ising spin glasses, i.e., that the leading 
correction-to-scaling exponent lo is large. More specifi- 
cally, in Ref. [B| a value for ui between 1.3 and 1.6 was 
found. Although the range of system sizes accessible to 
Monte Carlo simulations in four space dimensions is more 
limited than in three space dimensions, the model poses 
a "good compromise" case where corrections can be kept 
small while the system sizes are reasonably large. In this 
article we choose the observables that display the small- 
est corrections to scaling for the system sizes accessible 
via Monte Carlo simulations. We feel that introducing 
correction-to-scaling terms along the lines of Ref. [HI in 
order to give more precise estimates for the critical ex- 
ponents cannot be controlled sufficiently well, and thus, 
we do not display an analysis using these methods. How- 
ever, we show that our results are perfectly compatible 
with a large correction-to-scaling exponent uj. 

As with the three-dimensional Ising spin glass , 3 ' 16 ! 17 
there have been many different estimates for the critical 
exponents (especially for the anomalous dimension 77) of 
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the four-dimensional Ising spin glass as shown in Table 

HI 

The main conclusion of this work is that equilibrium 
universality in four-dimensional spin glasses is satisfied, 
since we find agreement within error bars for all the finite- 
size scaling functions and critical exponents studied. 

In Sec. |TI]we introduce the model as well as the mea- 
sured observables. In addition, we describe the numerical 
methods used. Results are presented in Sec. IIII) followed 
by concluding remarks in Sec. HVl Finally, a discussion 
of other commonly-used observables that have been less 
useful in the present analysis is given in Appendix[X] De- 
tails of the analysis are presented in Appendices [El and [C| 
as well as an extended scaling analysis of the data with 
Gaussian disorder in Appendix [Dl 



II. MODEL, OBSERVABLES, AND 
NUMERICAL DETAILS 

A. Edwards- Anderson model 



is dimensionlcss, and scales as 

g(L, T) = g[BL 1/u {T - T c )\ + corrections. (5) 

Here T c is the critical temperature and B is a metric fac- 
tor. The critical exponent v describes the divergence of 
the infinite- volume correlation length £(T) as the tem- 
perature approaches T c , i.e., £(T) ~ |T - T c \~ v . The 
corrections to scaling in Eq. ([5|), as well as in Eqs. ([9|) 
and (fi"2"|) are asymptotically dominated by the leading 
correction-to-scaling exponent u> and vanish in the ther- 
modynamic limit (L — > oo). In Eq. ((4]) (O) represents a 
thermal average of an observable O, [0] av is a disorder 
average, and q is the spin overlap; 
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In Eq. ([6]) "a" and "6" represent two replicas of the sys- 
tem with the same disorder. In addition, we study the 
finite-size correlation length ; 32 i 33 i 34 



The Edwards-Anderson (EA) Ising spin-glass 
Hamiltoniani^ is given by 



1 



2sin(|k min |/2) 



Xsg(O) 

XSG(kmin) 



1/2 



(7) 



H — ^2 JijSiSj 



(1) 



The Ising spins S, = ±1 lie on a hypcrcubic lattice of 
size N — L d in d = 4 space dimensions with periodic 
boundary conditions. The sum is over nearest neighbors 
on the lattice. We study two versions of the model: 

(i) Gaussian-distributed interactions Jy with zero 
mean and standard deviation unity; 
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2tt 



(2) 



(ii) Bimodal-distributed random interactions^ with a 
bond dilution of 65%, i.e., 

V(Jij)= (l-p)S(J ij ) 

+ f + 6(^ + 1)], 



with p = 0.35 (Rcf. [3 



where k m ; n = (2n/L, 0, 0, 0) is the smallest nonzero wave 
vector and the wave-vector-dependent spin-glass suscep- 
tibility is given by 



XsgI 



(k) = ^E[^^) 2 ] a v^ k - (Ri - Rj) - 

i,3 



(8) 



The finite-size correlation length £(L,T) divided by the 
system size is also a dimensionlcss quantity, i.e., 



corrections. 



(9) 



For L — > oo data for £(L,T)/L as well as for g(L,T) 
intersect at T — > T c . For finite systems the data cross 
at an effective critical temperature T* (L) that converges 
asymptotically to T c as^i 



T*{L)-T C oc L~ u ~ l l v . 



(10) 



In the following we denote the value of an observable O 
measured at this effective critical point by 



0*=0[T* c {L)}. 



(11) 



B. Measured observables 

In order to compute the critical parameters and hence- 
forth test for universality, we compute different observ- 
ables that are known to show a good signature of the 
phase transition. The Binder cumulant^i defined via 



?(L ' T) -2^ 3 wm 



(4) 



The definition of the finite-size correlation length £(X, T) 
in Eq. ([7]) involves in general the same leading correc- 
tions to scaling as xsGj which in turn is given by lu. 
Furthermore, this definition is not unique^ and differ- 
ent definitions of £(L,T)/L show differences of the order 
of L~ 2 . Such differences may seem irrelevant, but for 
the small systems sizes that can be accessed in spin-glass 
simulations correction terms might actually be visible in 
the data. Note that in Eqs. ©, ©, and ([12}. T c and 
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TABLE I: Different estimates (sorted chronologically and with respect to the method used for the determination) of the 
critical exponents computed by different groups for Gaussian (G), bimodal (±J), triangular (T), uniform (U) and Laplacian 
(L) random bonds. The estimates, especially the ones for 77, show strong variations and often do not agree. The last two rows 
show the results from this study. The critical temperatures denoted by an asterisk indicate that the variance of the coupling 
distribution used in the corresponding study is not normalized to unity (Ref. fish. 



Authors 


Couplings 


Method 


7 


V 


n 


T c 


Singh & Chakravarty (Ref. 19) 


±J 


Series 


2.0(4) 






2.02(6) 


Klein et al. (Ref. 20) 


±J 




2.00(25) 


0.95 


-0.11 


2.04(3) 


Daboul et al. (Ref. J3) 


U 




2.4(1) 






1.10(2)* 




±J 




2.5(3) 






1.96(7) 




G 




2.3(1) 






1.79(1) 




T 




2.5(2) 






1.36(3)* 


Bernardi & Campbell (Ref. _21) 


±J 


Dynamic MC 






-0.31(1) 


1.99(1) 




G 








-0.47(2) 


1.77(1) 




U 








-0.37(2) 


1.91(1) 




L 








-0.60(3) 


1.52(1) 


Bhatt & Young (Ref. .22) 


G 


Static MC 


1.8(4) 


0.8 


-0.30(15) 


1.75(5) 


Reger et al. (Ref. 23) 


±J 








-0.5 




Parisi et al. (Ref. 24) 


G 




2.1(2) 


0.9(1) 


-0.35(5) 


1.80(1) 


Ney-Nifle (Ref. 25) 


G 






0.87(15) 




1.80(5) 


Marinari & Zuliani (Ref. 26) 


±J 






1.0(1) 


-0.30(5) 


2.03(3) 


Hukushima (Ref. 27) 


±J 






0.92(6) 




2.00(4) 


This study 


G 


Static MC 


2.32(8) 


1.02(2) 


-0.275(25) 


1.805(10) 




±J bond-diluted 




2.33(6) 


1.025(15) 


-0.275(25) 


1.0385(25)* 



the metric factor B are nonuniversal, but, since B is in- 
cluded explicitly, the scaling functions g(x) and t;(x) are 
both universali^^ Since for both disorder distributions 
studied we use the same boundary conditions and sam- 
ple shapes, these scaling functions are expected to be the 
same for different disorder distributions if the systems are 
in the same universality class. This is a necessary yet not 
sufficient condition. In addition, the critical exponents, 
as well as the values of the scaling functions at criticality 
(T = T c ) have to agree. 

In addition to studying the Binder cumulant as well 
as the finite-size correlation length (from which we ob- 
tain T c and the critical exponent v), we need to study 
another observable to obtain a second critical exponent 
to fully characterize the universality class of the models 
Thus we also study the scaling behavior of the spin-glass 
susceptibility xsg = Xsc(k = 0) (also xsg = N[(q 2 )] av ). 
The spin-glass susceptibility scales as 

Xsg{L,T) = CL 2 -"x[ J BL 1 / !/ (T-T c )]+corrections, (12) 

where the anomalous dimension rj is the second critical 
exponent needed to establish the universality class of the 
model. In Eq. Q12[) C represents a nonuniversal ampli- 
tude. We also study another related quantity £(L,T) 
that has shown to be useful in determining the critical 
exponent 77, which is defined as2£ 



C(L,T) = ?p. 



(13) 



The advantage of studying £ is mainly given by the fact 
that the statistical correlations between xsg an d £, lead 
to smaller errors on £ than on xsg an d therefore to a 
more precise determination of rj. 



Finally, we study other phenomenological couplings 
that have been suggested to compute the spin-glass tran- 
sition temperature, which are the lack of self-averaging 
A given by^MOdid^ 



A(L, T) = K« 2 ) 2 l av 



K<7 2 )] 



2\12 
av 



[<<? 



2\12 
av 



and the Guerra parameter G given by*« 



44.45 



G{L,T) 



[(9 



2\2l 



- K<? 2 )] 2 



[<<7 4 >]av-[(<Z 2 > 



2\12 



(14) 



(15) 



Both A and G are related to the Binder cumulant g 
through the relation g = 1 — A/ (2G). 

The scaling expressions in Eqs. (0, @, and (TT2"j) can 
be used to determine the critical exponents, but in prac- 
tice this strategy is not very promising because especially 
in the case of the spin-glass susceptibility analytic correc- 
tions to scaling are very easily confused with the leading 
scaling behavior, which in turn leads to unreliable esti- 
mates of the critical exponents!^ In order to determine 
reliable estimates for the critical quantities, we use the 
quotient method, which avoids the problems of the ana- 
lytic corrections to scaling in an elegant manner i 34 ' 46 i 47 
For any observable 0(L,T) and the finite-size correlation 
length, finite-size scaling theory predict a 48 ' 49 i 50 that 



Q(L,T) 
C(oo, T) 

as well as 

Q(sL,T) 
0(L,T) 



o 



tt°°,T)/L 



&L,T)/L; 



(16) 



O 



(r w ,zr"), (17) 
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where fo and Fq are universal finite-size scaling func- 
tions and s > 1 is a scale factor. The exponent uj is 
again the leading nonanalytic correction-to-scaling expo- 
nent. Because in Eq. (jf 7[) only pairs of finite system sizes 
L and sL appear, this formulation is well adapted for use 
in numerical simulations. For example, the knowledge 
of the universal scaling functions F x and F% [meaning 
O = x or O = £ in Eq. (fTT|) . respectively] allows us to 
extract the critical exponents r/ and v using the quotient 
method. For the quotient method one defines an effective 
critical temperature T* at which the correlation length 
measured in units of the lattice size L is equal for the 
pair of systems, i.e., 



£(L,T*)/L = T*)/(sL), 



or alternatively, 



g{L,T* c )=g{sL,T* c ). 



(18) 



(19) 



Note that we do not use a different notation for T* de- 
fined through Eqs. (p~8|) and (p~9|) although the corre- 
sponding T* in general is different, and only in the ther- 
modynamic limit converges to a unique value T c as indi- 
cated in Eq. (flTj|) . The crossings are determined by fitting 
a cubic spline through the data. We have refrained from 
using rewcighting techniques for the determination of the 
crossings ) 46 i 51 because the remaining statistical errors of 
the sample average dominate the errors even for the good 
statistics we have at hand. In addition to Eq. (|17p we also 
study the following equivalent relation 



Q{sL,T) 
G(L,T) 



g{L,T);s 



(20) 



which, however, will prove to be advantageous in the 
present study. Another case in which this version of the 
finite-size scaling relation is very useful is to study scaling 
properties within the spin-glass phased 

The critical exponent xq associated with a given ob- 
servable which at criticality diverges as (T — T c )~ x ° can 
then be estimated via the quotient 



jo/" 



Q(sL,T*) 
0(L,T*) 



+ 0(L~"), 



(21) 



and thus from the finite-size scaling function Fq. In 
Eq. (|21[) the critical exponent v of the correlation length 
is unknown, but it can be estimated for example from the 
finite-size scaling function of the temperature derivative 
of the correlation length £(L,T), Fg T £, via 



V» = l + ^d x Fz(x,s)\+0(L-"), 



(22) 



with x = £(L,T)/L and x* = £(L,T*)/L. Alternatively, 
v can also be estimated from the temperature derivative 
of the finite-size scaling function of the Binder cumulant 

FdTgi 1*6., 

s 1 '" = 1 + g*d g F g (g, s)\ g=g , + O(L-). (23) 



Here g* = g(T*,L). In our study we fit cubic splines 
to the data of F% and F g to calculate the derivatives in 
Eqs. ([22]) and (p3]l. A detailed derivation of Eqs. |22|) 
and (|23[) is given in Appendix IB"] 

The anomalous dimension rj can be determined from 
F x . Using the scaling relation 



we obtain 



7= i/(2 - rj), 



s 2 -" = F* + 0{L~ U ). 



(24) 



(25) 



Similarly, rj can be obtained from the finite-size scaling 



function F, 



-') _ 



(26) 



We determine the critical exponent 7 from Eq. (|24[) using 
the estimates of v and r\. Finally, we also study the 
Binder cumulant as a function of £(L,T)/L (Ref. 0i, 
and I53T ) . Nonuniversal metric factors cancel and so 
g(L,T) = g[£(L,T)/L] with g a universal function.—^ 
Therefore data for different disorder distributions should 
fall on the same universal curve if the models share the 
same universality class. 



C. Simulation details 

The simulations are done using exchange (parallel tem- 
pering) Monte Carl o 54 ' 55 and the simulation parame- 
ters are presented in Tables [XT] and IIIII for the Gaussian 
and bond-diluted bimodal disorder distributions, respec- 
tively. For the Gaussian disorder we test equilibration 
using the method presented in Ref. [H. For the bond- 
diluted bimodal disorder we use a multispin coded ap- 
proach to speed up the simulation in addition to a clus- 
ter updating routine^i 5 ^ which can substantially speed 
up equilibration when the system is diluted and com- 
plements the parallel tempering Monte Carlo updates. 
Although we did no systematic study, we have the im- 
pression that in four dimensions the additional cluster 
updates are somewhat less effective than in lower space 
dimensions. Equilibration in the bond-diluted case is 
tested by a logarithmic binning of the data. Once the 
last three bins for all obscrvables agree within error bars, 
we declare the system to be in thermal equilibrium. 



III. RESULTS 

In the following we perform a finite-size scaling anal- 
ysis of different obscrvables for both models. We deter- 
mine quantities with small scaling corrections that pro- 
vide good estimates for the critical temperature and criti- 
cal exponents. We then compare different finite-size scal- 
ing functions for the model with different disorder distri- 
butions and discuss the influence of corrections to scaling 
on our results. 
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T T 

FIG. 1: (Color online) Determination of the critical temperature of the four-dimensional Edwards- Anderson Ising spin glass 
with Gaussian disorder. In panel (a) the finite-size correlation length £(L,T)/L as a function of the temperature T for different 
system sizes L is shown. The data for L > 6 cross at T c = 1.810(15). In panel (b) the corresponding data for the Binder 
cumulant g(L,T) as a function of the temperature T for different system sizes L are shown. The data cross at T c — 1.805(10), 
in agreement with the data for the correlation length. The crossing of the data for the Binder cumulant is cleaner than for the 
correlation length and shows no noticeable drift of the crossing point with increasing system sizes. 
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FIG. 2: (Color online) Determination of the critical temperature of the four-dimensional Edwards-Anderson Ising spin glass 
with diluted bimodal disorder. In panel (a) the finite-size correlation length £(L,T)/L as a function of the temperature T for 
different system sizes L is shown. The data for L > 6 cross at T c = 1.042(5). In panel (b) the corresponding data for the 
Binder cumulant g(L, T) as a function of the temperature T for different system sizes L are shown. The data cross cleanly at 
T c = 1.0385(25). The crossing of the data for the Binder cumulant is again more precise than the one of the correlation length 
and shows no noticeable drift of the crossing point with increasing system sizes. 



A. Gaussian disorder (no dilution) 

Panel (a) of Fig. [T] shows the data for the finite-size 
correlation length as a function of the temperature for 
different system sizes L. With increasing L, the data 
show a shift of the effective T c toward a smaller value of 
T c . This effect has already been observed in studies of the 



three-dimensional EA model,— however, in contrast to 
the three-dimensional case the range of available lattice 
sizes in four dimensions is more restricted and therefore 
the effect of this shift is clearly a restriction for a pre- 
cise determination of T c using the crossings of £(L,T)/L. 
Taking the data for L > 6 and neglecting the remaining 
scaling corrections, we obtain T c = 1.810(15). 



6 



2.4 
2.2 



Ex 



2 — 



1.4 
1.2 
1 



i i | i i i | i i 
(a) bond-diluted bimodal © * 



O 6/3 - 
A 8/4 
□ 10/5 



0.2 



0.4 



0.6 



f(L.r)/L 




FIG. 3: (Color online) Finite-size scaling functions of the correlation length F% and the Binder cumulant F g of the four- 
dimensional Edwards- Anderson Ising spin glass with bond-diluted bimodal disorder for different system sizes L. The comparison 
between the scaling function F% as a function of the finite-size correlation length f(L, T)/L is shown in panel (a) and F g as a 
function of the Binder cumulant g(L, T) shown in panel (b). Given the small system sizes studied, the Binder cumulant is more 
suited as a scaling variable, since it displays clearly smaller scaling corrections than the finite-size correlation length. The data 
collapse for the scaling function F g is excellent for the small system sizes studied. The broken lines in both panels indicate the 
condition that defines the effective critical point T*. 



TABLE II: Parameters of the simulations for the Gaussian- 
distributed disorder. L denotes the system size. N sa , is the 
number of samples and iV sw is the total number of Monte 
Carlo sweeps performed on a single sample for each of the 
2Nt replicas. T m i n and T max are the lowest and highest tem- 
peratures simulated, and Nt is the number of temperatures 
used in the parallel tempering method. 



L 


N sa 


iVsw 


N T 




Tmax 


3 


20000 


131072 


29 


1.400 


3.061 


4 


20000 


131072 


29 


1.400 


3.061 


5 


20000 


131072 


29 


1.400 


3.061 


6 


20000 


131072 


29 


1.400 


3.061 


8 


3500 


524288 


29 


1.400 


3.061 


10 


2000 


524288 


29 


1.400 


3.061 



TABLE III: Parameters of the simulations for the bond- 
diluted bimodal disorder distribution. For details see Table 



mi 



L 


iVsa 


iVsw 


N T 


-Maxin 


Tmax 


3 


102400 


20000 


11 


0.950 


1.800 


i 


107680 


40000 


11 


0.950 


1.800 


5 


101699 


40000 


11 


0.950 


1.800 


6 


101664 


40000 


11 


0.950 


1.800 


8 


41408 


100000 


21 


0.950 


1.800 


10 


24160 


100000 


21 


0.950 


1.800 



Panel (b) of Fig. Q] shows the data for the Binder cumu- 
lant g(L, T) as a function of the temperature for different 
system sizes. The data cross cleanly at the critical tem- 
perature of T c — 1.805(10) and there is no shift as for 



the crossings of the finite-size correlation length [panel 
(a)]. In contrast to the situation in three dimensions the 
data splay out very well below T c making a precise de- 
termination of T c from the Binder cumulant data possi- 
ble, although the error bars on the Binder cumulant are 
slightly larger than the ones on the correlation length. 
We consider Fig. [T] as a first indication that it might be 
profitable to use g(L,T) instead of £(L,T)/L as scaling 
variable in the finite-size scaling analysis in four dimen- 
sions, possibly because T c is considerably larger and thus 
the crossing point further away from T = where g — > 1. 

In Table HVl ( Appendix ICj) we present the results for the 
critical quantities we obtain from the quotient method 
using T* defined from the crossings of £/L and g, respec- 
tively. While the results for T c , g(T c ), and v show no 
noticeable scaling corrections, the ones for 77, £(L,T C )/ L, 
and in a minor extent 7 indicate clearly the presence of 
such corrections. In conclusion we obtain the following 
values for the critical quantities for the undiluted Gaus- 
sian disorder: 

T c = 1.805(10), 
g(T c ) - 0.470(5), 
v = 1.02(2), 
77 = -0.275(25), 

7 = 2.32(8). (27) 

Our value for 77 contains a crude extrapolation to the 
thermodynamic limit, which makes use of the fact that 
corrections to scaling seem to be disappearing very fast 
with increasing system size. The result we give for i] 
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open symbols bond-diluted bimodal 4? 
0.6 — full symbols Gaussian w 




0.2 0.4 0.6 



f(L,T)/L 

FIG. 4: (Color online) Comparison of the Binder cumulant g 
as a function of the finite-size correlation length for the four- 
dimensional EA Ising spin glass with Gaussian (full symbols) 
and diluted bimodal disorder (open symbols) for different sys- 
tem sizes L > 8. Both functions agree very well; further evi- 
dence for universal behavior. 

is justified later in Fig. [5J The value for 7 is deter- 
mined from our estimates of v and 77. Its value depends 
much more on a precise determination of v than on a 
precise determination of ?y. We did not try to deter- 
mine an infinite-volume extrapolation for the value of 
£(L,T C )/L, which is another universal quantity, because 
the data from such a limited range of system sizes do 
not allow for a controlled extrapolation. Note that in 
the case of 77, where one might expect similar problems, 
we rely on the convergence of the estimates from oppos- 
ing sides for the different scaling functions (see Fig. 
thus allowing for a somewhat more reliable extrapola- 
tion. Taking into account the presence of corrections to 
scaling, the estimates for the critical exponents obtained 
from the extended scaling method of Campbell et al^ 
(see Appendix [D)) do agree with our final estimates. 



ing than the correlation length. This is confirmed by a 
comparative analysis of the finite-size scaling functions 
Fg as a function of £(L,T)/L and F g as a function of 

g(L,T). 

The finite-size scaling functions Fg and F g contain in 
principle the very same information. However, due to the 
fact that we are working in a restricted range of lattice 
sizes there are visible differences. In panel (a) of Fig. [3] 
we show the scaling function Fg and in panel (b) of Fig. [3] 
of the corresponding scaling function F g . Comparing the 
two figures one clearly sees that F g shows much smaller 
finite-size corrections than Fg. This figure is clear evi- 
dence that for our data the Binder cumulant g(L,T) is 
better suited as a scaling variable than the correlation 
length £(L,T)/L. This is in contrast to recent work in 
three dimensions, where typically £(L,T)/L is used in- 
stead of g(L,T) in order to cleanly determine T c and 
the critical exponents^^ The reason for this behavior is 
possibly given by the following arguments: First, in three 
dimensions the crossing of the Binder cumulant data is 
at a rather flat angle making it difficult to determine T c 
precisely, while in four dimensions the crossing is at a 
much steeper angle. Second, the range of available sys- 
tem sizes in four dimensions is smaller than in three di- 
mensions and therefore the use of observables with small 
scaling violations is better in order to have good control 
over the estimates of the critical quantities. 

In Table IIVI (Appendix [C]) we present the results for 
the critical quantities we have obtained from the quotient 
method using T* , defined either from the crossings of £/L 
or g, respectively. We obtain the following values for the 
critical quantities for the link-diluted bimodal disorder 
distribution: 

T c = 1.0385(25), 
g(T c ) = 0.472(2), 
v = 1.025(15), 
7/ = -0.275(25), 

7 = 2.33(6). (28) 
The estimate of 77 is justified below in Fig. [5J 



B. Diluted bimodal disorder 

Panel (a) in Fig. [2] shows data for the finite-size cor- 
relation length as a function of the temperature for dif- 
ferent system sizes. With increasing system size L the 
data show — as for the Gaussian data — a noticeable shift 
of the effective critical temperature T c toward a smaller 
value. Taking the data for L > 6 and neglecting the 
remaining scaling corrections we obtain T c = 1.042(5). 
Panel (b) in Fig. [2] shows the data for the Binder cumu- 
lant g(L, T) as a function of the temperature for different 
system sizes. The data cross cleanly, and we determine 
T c = 1.0385(25). We find again as in the Gaussian case 
that the Binder cumulant has smaller corrections to scal- 



C. Comparison of finite-size scaling functions 

The results for the critical quantities given in Eqs. ([27]) 
and (|28|) are consistent with universal critical behavior of 
the four-dimensional Edwards- Anderson Ising spin-glass 
model. We further strengthen this result by comparing 
the finite-size scaling functions of the two models not 
only at the critical point but within the whole scaling 
region. The direct comparison of the finite-size scaling 
functions has shown to be probably the best approach 
to check for universality in spin glasses as it allows for a 
completely parameter-free comparison of different mod- 
els (see Refs. [H and[H for a comparison of different mod- 
els in three space dimensions). In Fig. 3] we compare 
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FIG. 5: (Color online) Comparison of the finite-size scaling function F g shown in panel (a) and F x shown in panel (b) for 
the four- dimensional model with Gaussian and bond-diluted bimodal disorder for different system sizes I as a function of the 
Binder cumulant g. The data for both finite-size scaling functions show very little scaling corrections. The fact that the curves 
for the two different models for F g and F x fall on one single curve is a strong evidence for universal critical behavior of the 
four-dimensional Edwards- Anderson model. The insets in both panels present an enlarged view around the critical point [which 
is at g(L,T c ) = 0.472(2)] and show in more detail that the data for the finite-size scaling functions collapse onto one single 
master curve. The data point with the label "oo" in panel (b) indicates our infinite-volume extrapolation of F x at criticality 
corresponding to g = 0.275(25) and g(T c ) = 0.472(2). 
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spectively, at T* (defined by £,* / L and g*) as a function of 
with uj = 2.5. The data for the Gaussian (full symbols) 
and bond-diluted bimodal (open symbols) disorder are consis- 
tent within errors and they are also consistent with a unique 
value of rj = —0.275(25) for an infinite-volume extrapolation. 



the Binder cumulant g(L,T) plotted against the correla- 
tion length £(L,T)/L for the two different disorder dis- 
tributions. The data agree and collapse onto a single 



master curve, which is a strong evidence for universality. 
In panel (a) of Fig. [5] we compare the finite-size scaling 
function F g as a function of the Binder cumulant g(L, T). 
The data collapse again within error bars onto a single 
master curve. This result illustrates once more that the 
estimates of the critical exponent v of the two disorder 
distributions do agree within error bars. In panel (b) of 
Fig. [5] we compare the finite-size scaling function F x as a 
function of the Binder cumulant g{L,T). Again the data 
collapse is within the error bars, which means that also 
the estimates of the critical exponent rj of the two models 
coincide within errors. 

Finally, in Fig. [6] we show the convergence of the ef- 
fective exponent T] e s from F x and through Eqs. 
and (j26|) . respectively, measured at the crossings of £/L 
and g as a function of All the scaling corrections 

are compatible with a leading correction-to-scaling expo- 
nent of u> « 2.5. Clearly, this is only an effective ex- 
ponent because for the small system sizes we have at 
hand, we cannot expect to be in the asymptotic scaling 
regime where the leading correction-to-scaling exponent 
dominates. However, our estimate of uj is consistent with 
the assumption that the corrections to scaling we see are 
dominated by the first nonleading correction-to-scaling 
term since our value for uj is roughly twice as large as the 
one obtained from high-temperature expansion studies 
in Ref. or alternatively they might be due to remain- 
ing analytic corrections with an effective correction-to- 
scaling exponent given by 2 — rj. Figure [5] shows that an 
infinite- volume extrapolation of the different estimates is 
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FIG. 7: (Color online) The Guerra parameter G(L, T) [panel (a)] and the lack of self- averaging parameter A(L, T) [panel 
(b)] as a function of temperature T for the four- dimensional Edwards- Anderson model with bond-diluted bimodal disorder for 
different system sizes L. The crossings shift noticeably to a smaller effective T c in both cases for increasing system size making 
an accurate determination of T c from these data difficult and imprecise. In the case of A(L, T) the fact that the crossing occurs 
close to the maximum of the curves makes the situation for an accurate determination of T c (vertical dashed lines) impossible. 
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FIG. 8: (Color online) The Guerra parameter G(L, T) [panel (a)] and the lack of self-averaging parameter A(L, T) [panel (b)] 
as a function of the Binder cumulant g(L, T) for the four-dimensional Edwards- Anderson Ising spin glass with bond-diluted 
bimodal disorder for different system sizes L. Both G(L, T) and A(L, T) display large corrections to scaling. 



compatible with a unique value of 77 = —0.275(25) in the 
infinite-volume limit for both models. 



IV. CONCLUSIONS 

We have tested universality in four-dimensional Ising 
spin glasses and computed precise estimates of the crit- 
ical parameters of the model with both Gaussian and 
link-diluted bimodal disorder. Our results show that the 
different critical exponents for different disorder distri- 



butions agree well. Furthermore, by plotting the Binder 
ratio as a function of the correlation length, we show 
that four-dimensional Ising spin glasses (with compact 
disorder distributions) seem to share the same universal- 
ity class. Furthermore, we compute different finite-size 
scaling functions in four space dimensions defined via ra- 
tios of different observables and show that these show 
small corrections to scaling, especially when studied as 
a function of the Binder parameter. The results pre- 
sented thus indicate that universality is not violated in 
four-dimensional spin glasses. 



10 



10' 



10° 



10- 



x 



io- 



io- 



Gauss 



° 9AB9, 



L 

O 3 
A 4 

' □ 5 
O 6 
O 8 
OlO 




1.805 
rf = -0.275 
v = 1.03 



10-' 



10 8 



10° 10' 

(LT) l /"[l - (T c /T) 2 ] 

FIG. 9: (Color online) Scaling plot of the spin-glass suscep- 
tibility according to the extended scaling approach for the 
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APPENDIX A: OTHER OBSERVABLES 

In this appendix we discuss the quantities that have 
been shown to be less useful in the study of the lo- 
cation of the spin-glass transition and the issue of the 
universality of the four-dimensional Edwards-Anderson 
model. We have decided to present our results concern- 
ing these quantities since our data for the bond-diluted 
bimodal coupling distribution has by far the best statis- 
tics in the context of Monte Carlo simulations of the four- 
dimensional Edwards- Anderson model. 

Panel (a) of Fig. [7] shows the data for the Guerra pa- 
rameter G(L, T) defined in Eq. (|15|) as a function of 
the temperature for different system sizes. The cross- 



ings of the data for increasing system sizes shift notice- 
ably toward smaller temperatures. We find that G(L, T) 
has rather large finite-size corrections and also relatively 
large errors compared to, e.g., the Binder cumulant and 
therefore is not well suited for an accurate determination 
of T c . The same conclusions have been found by Balles- 
teros et al£^ for the three-dimensional bimodal Ising spin 
glass. Note, however, that in the case of mean-field spin 
glasses the situation may be different as the Guerra pa- 
rameter has shown to be more efficient in locating a 
spin-glass transition than the Binder cumulant in certain 
situations^ 

Panel (b) of Fig. [7] shows the data for the lack of self- 
averaging parameter A(L, T) defined in Eq. (fl"4|) as a 
function of the temperature for different system sizes. 
This parameter is related to the Guerra parameter and 
also shows large finite-size corrections and due to the fact 
that the crossings happen close to the maximum of the 
curves, where the slope changes very fast, the crossings 
(apart from the fact that they move noticeably) cannot 
be determined reliably. This fact together with the rather 
large relative error makes that this quantity the least 
suited for a precise determination of the critical temper- 
ature of all the quantities discussed here. 

In Fig. [8] we show the data for G(L,T) [panel (a)] and 
A(L,T) [panel (b)] as a function of the Binder cumu- 
lant. Clearly, the two quantities have strong finite-size 
corrections. The fact that these quantities may present 
such strong finite-size scaling corrections should be kept 
in mind when, e.g., the behavior of A{L,T) is used to 
make statements on the nature of the spin-glass phase as 
it is done in Ref. l58l. 



APPENDIX B: DERIVATION OF THE 
QUOTIENT RELATIONS 



In the following, we derive Eqs. ([22]) and (|23[) in detail. 
Starting from the definition of we can write, 



asL,T) = F^(L,T)/L;s]aL,T). 



(Bl) 



Taking the derivative with respect to T, using the chain 
rule and finally dividing by c?t£(L, T), we arrive at 



&r£(sL,T) _ £(L,T) 
dr£{L, T) L 

Note that 



d x F c (x;s)+F^(L,T)/L]. (B2) 



F aT6 [S(L,T)/L;s] = 



d T gsL,T) 



(B3) 



Using the fact that at the effective critical point, where 
£(£, T*)/L = x* , we have F^(x*; s) = s and the fact that 
the correlation length £ close to the critical point has a 
simple scaling form given in Eq. ([9]), we see that 



F dTi (i(L,T)/L- S ) 



■0(L U ) 



x*d x F 6 (x;a)\ iB=xm +a, (B4) 
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TABLE IV: Results from the quotient method for the Gaussian (top) and link-diluted bimodal (bottom) disorder at T* 
determined from the crossings of £/L and g, respectively. We group the data according to the observables which were used to 
compute them (separated by vertical bars). 



Crossings of f/L L 


T* 




F* r]cS 7eff 


F£ rj cff 7 eff 


3/6 
4/8 
5/10 


1.858(6) 
1.824(7) 
1.814(8) 


0.396(2) 5.00(3) 1.007(8) 
0.418(3) 4.67(5) 1.018(14) 
0.427(4) 4.62(6) 1.010(17) 


4.52(6) -0.18(2) 2.20(3) 
4.68(9) -0.23(3) 2.27(6) 
4.74(12) -0.24(4) 2.26(8) 


1.129(7) -0.175(9) 2.19(3) 
1.168(11) -0.224(16) 2.26(5) 
1.184(12) -0.244(15) 2.27(5) 


Crossings of g L 


T* 


9* d g F* u eS 


F* 77eff 7eff 


F£ rj cff 7eff 


3/6 
4/8 
5/10 


1.803(6) 
1.805(8) 
1.805(8) 


0.470(3) 2.01(5) 1.042(25) 
0.469(5) 2.08(9) 1.018(42) 
0.471(7) 2.06(11) 1.022(54) 


5.14(7) -0.36(2) 2.46(8) 
4.96(11) -0.31(3) 2.35(13) 
4.92(14) -0.30(4) 2.35(17) 


1.160(9) -0.214(12) 2.31(7) 
1.186(14) -0.246(17) 2.29(11) 
1.198(20) -0.261(24) 2.31(15) 




Crossings of £/L L 


' c 




F* r]cS 7eff 


F£ r/ c fi 7eff 


3/6 
4/8 
5/10 


1.0716(15) 
1.0502(16) 
1.0441(18) 


0.400(2) 4.86(2) 1.021(7) 
0.424(2) 4.56(2) 1.025(7) 
0.433(2) 4.42(3) 1.032(9) 


4.51(4) -0.17(1) 2.22(3) 
4.68(6) -0.23(2) 2.29(3) 
4.74(7) -0.25(2) 2.32(4) 


1.127(6) -0.172(8) 2.22(3) 
1.170(6) -0.227(8) 2.28(3) 
1.186(6) -0.246(8) 2.32(3) 


Crossings of g L 


T* 
J- c 


g* d g F* u c g 


F* ri c g 7 e ff 


F£ r/ c g 7 e ff 


3/6 
4/8 
5/10 


1.0390(18) 
1.0390(18) 
1.0384(20) 


0.472(2) 2.09(3) 1.010(14) 
0.472(2) 2.06(3) 1.021(15) 
0.473(3) 2.06(4) 1.019(20) 


5.12(4) -0.36(1) 2.38(3) 
4.94(5) -0.31(2) 2.36(5) 
4.92(6) -0.30(2) 2.34(7) 


1.163(5) -0.224(7) 2.25(3) 
1.185(5) -0.245(7) 2.29(4) 
1.197(6) -0.259(8) 2.30(6) 



from which we obtain Eq. (|22|) by dividing by s. The 
derivation of Eq. (|23|) has as the starting point 

g(sL,T) = F g (g(L,T);s)g(L,T), (B5) 

and for the rest is analogous to the one given for Eq. (f!?2]) . 

APPENDIX C: RESULTS FROM THE 
QUOTIENT METHOD 

In this section wc list the detailed results from the quo- 
tient method. The data are grouped by the observable 
used to compute the estimates in Table HVl 

APPENDIX D: EXTENDED SCALING 

Recently, Campbell et al^ suggested an extended 
scaling approach, which allows one to extend the scaling 



region from \L^ V (T -T c )\ < 1 to virtually T -> oo. The 
method has the advantage in that it drastically reduces 
the corrections to scaling commonly found when perform- 
ing a simple finite-size scaling analysis of the spin-glass 
susceptibility— In that scaling approach the scaling equa- 
tion for the susceptibility [see Eq. ([12]) ] is modified in the 
following ways^ 



XsG (i,T) ~ (LT) 2 ~ V x[B(LT)^\l - {T c /T)% (Dl) 



In Fig. [9] we illustrate the quality of the critical param- 
eters by performing an extended hnite-size scaling plot 
of the susceptibility for the EA spin glass with Gaussian 
disorder obtained from the quotient method i 34 ! 46 Simi- 
lar results are obtained for the model with bond-diluted 
bimodal disorder, as well as other observables. 
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